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ABSTRACT 

N 

Analytical  and  computer  approaches  are  presented  for  predicting  the  site 
errors  of  airborne  crossed  loops  used  for  lightning  direction-finding.  This 
type  of  azimuth  error  is  caused  by  distortion  of  the  incident  field  due  to 
induced  skin  currents.  Computer  code  results  are  presented  for  a T-39  and  a 
C-130  aircraft  used  in  f light-testing  a commercial  cb£paed-loop  system.  Another 
type  of  bearing  error,  caused  by  non-vertical  lightning  channels,  is  analyzed 
for  the  case  where  both  detector  and  lightning  source  are  above  a conducting 
ground  plane.  Also  discussed  is  the  180-degree  ambiguity  and  its  resolution 
by  an  electric  field  antenna. 


AIRBORNE  CROSSED  LOOPS  may  be  used  to  sense  the 
azimuth  of  lightning  discharges  as  an  on-board 
aid  in  the  avoidance  of  thunderstorm  regions. 
However,  the  loops  can  be  susceptible  to  serious 
site  errors  in  lightning  direction-finding, 
depending  on  the  geometry  of  the  aircraft  and 
the  loop  location  on  it  (1,  2,  3)*.  These  site 
errors  are  caused  by  scattering  (re-radiation 
due  to  induced  eddy  currents)  of  the  incident 
magnetic  wave  field  by  nearby  conducting  sur- 
faces (the  skin  of  the  aircraft).  These  eddy 
currents  produce  a secondary  magnetic  field  that 
results  in  a total  magnetic  vector  of  generally 
different  (distorted)  orientation  from  the  inci- 
dent orientation.  The  crossed-loop  detector  is 
"misdirected"  when  it  senses  the  distorted 
orientation  and  indicates  an  erroneous  bearing. 
(The  term  "misdirection"  was  used  earlier  in 
this  sense,  and  may  have  been  coined  by  A.  Som- 
morfeld  (13).)  Another  consequence  of  the  site 
error  is  that  it  disturbs  also  lightning  ranging 
systems  that  use  relations  between  magnetic 
amplitude  and  distance,  since  the  amplitude  as 
well  as  the  direction  is  distorted  by  the  scat- 
tering effect. 

Two  approaches,  analytical  and  numerical, 
are  presented  for  analyzing  and  predicting  the 
site  errors  of  airborne  crossed  loops.  In  both 
approaches  we  assume  that  the  aircraft  is  small 
compared  with  the  radiation  wavelength  (magneto- 
static limit),  and  solve  Laplace's  equation  to 
obtain  the  scattered  magnetic  field  vector,  for 
arbitrary  azimuth  angles  of  incidence  of  the 
wave  from  the  lightning  radiation  source.  The 
site  error  arises  from  the  fact  that  the  indi- 
vidual orthogonal  magnetic  field  components  are 
changed  by  different  factors.  The  error  anal- 
yses presented  here  and  in  Refs.  (1)  and  (2)  do 
not  appear  to  have  been  published  previously. 

In  the  analytic  approach,  fuselages,  wings 
and  tails  are  modeled  by  ellipsoids.  An  impor- 
tant analytical  result  is  that  large  errors  in 
lightning  direction-finding  can  occur  (e.g.,  of 
the  order  of  20°)  even  in  the  case  where  the 
instrument  is  symmetrically  located,  e.g.,  cen- 
tered on  a long  fuselage.  The  error  can  be 
still  larger  if  the  instrument  is  mounted  near 
the  nose  or  tail,  or  near  edges.  (See  Ref.  3.) 

In  the  numerical  approach,  the  aircraft 
geometry  is  modeled  realistically,  and  numeri- 
cal solutions  are  obtained  using  a 3-D  computer 
code.  Computational  results  are  presented  for 
two  actual  aircraft,  a T-39  Sabreliner  and  a 
C-130  Hercules.  These  aircraft  have  been  used 
in  flight  tests  to  evaluate  a coramercially- 
available  crossed-loop  system  for  severe 
weather  avoidance  (4,  5,  5).  The  predicted 
errors  (of  the  order  of  10°  for  the  T-39  and 
20°  for  the  C-130)  are  consistent  with  flight- 
test  results  for  the  particular  mounting  loca- 
tibns  chosen.  The  calculations  indicate  that 
other  locations  would  reduce  the  error. 

By  a combination  of  theory  and  experiment, 
one  may  determine  correction  factors  for  site 
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errors  affecting  crossed  loops,  for  any  given 
airplane  geometry.  The  results  of  this  deter- 
mination would  suggest  optimum  locations  for 
the  placement  of  the  loops  on  the  aircraft. 

The  correction  factors  need  be  determined  only 
once.  If  this  is  done,  the  site  errors  can  be 
completely  eliminated  by  suitable  adjustment  of 
the  electronic  amplification. 

Also  discussed  are  azimuth  errors  due  to 
non-vertically-polarized  lightning  channels, 
and  the  180-degree  ambiguity  inherent  in 
crossed  loops.  A non-vertical  (or  "slant") 
lightring  channel  will  cause  a crossed  loop  to 
determine  an  erroneous  direction,  if  the  loop 
and  lightning  channel  are  at  different  alti- 
tudes. The  effects  of  loop  altitude  above  a 
conducting  ground  plane  (as  opposed  to  light- 
ning altitude)  on  azimuth  (or  bearing)  errors 
do  not  appear  to  have  been  considered  in  the 
literature.  An  analytical  formula  is  derived, 
generalizing  results  obtained  for  zero  loop 
altitude  (9,  10)  to  the  case  of  finite  altitude. 

The  use  of  an  electric  antenna  in  conjunc- 
tion with  the  loop  to  resolve  the  180-degree 
ambiguity  by  correlating  the  phases  of  the  elec- 
tric and  magnetic  vectors  is  an  important  prob- 
lem, and  is  considered  in  the  Appendix.  It  is 
clear  that  low  signal-to-noise  ratios  or  unin- 
tentional phase  shifts  in  the  electronic  ampli- 
fiers can  produce  180-degree  errors. 

ANALYTICAL  APPROACH 

Evaluation  of  the  scattering  of  incident 
electromagnetic  waves  is  generally  a difficult 
problem,  even  for  perfectly  conducting  bodies 
of  simple  shape,  e.g.,  spheres  (7).  If  we 
assume  that  the  radiation  wavelength  is  larger 
than  the  dimension  of  the  airplane,  the  problem 
is  simplified  because  Maxwell's  equations  are 
replaced  by  the  Laplace  equation  of  magneto- 
statics. This  approximation  is  valid  for  fre- 
quencies below  the  first  aircraft  resonance, 
i.e.,  for  frequencies  below  about  one  MHz  (8). 

In  the  magnetostatic  limit,  the  boundary  condi- 
tion at  the  aircraft  surface  is  that  the  mag- 
netic field  be  tangent  to  the  surface  (zero 
normal  gradient,  a Neumann  boundary  condition). 
That  is,  the  magnetic  field  is  excluded  from 
the  interior  of  the  aircraft  by  the  induced 
skin  currents.  This  boundary  condition  is  jus- 
tifiable based  on  the  fact  that  the  skin  depth 
(given  by  66/v1*  in  mm  for  copper,  where  v is 
the  frequency  in  Hz  (7))  is  only  0.3  mm  at  the 
frequency  50  kHz  of  interest  here,  and  is  there- 
fore less  than  typical  aircraft  skin  thicknesses 
(at  least  one  mm). 

In  magnetostatics  the  magnetic  field  vec- 
tor at  any  point  may  be  expressed  as  the 
gradient  of  a scalar  potential  function  4, 
where  f satisfies  the  Laplace  equation  in  the 
region  exterior  to  the  airplane,  has  a vanish- 
ing normal  derivative  on  the  airplane  surface, 
and  asymptotically  approaches  a linear  form 
(constant  gradient)  at  infinity. 

For  analytical  purposes  (as  opposed  to  the 
numerical  approach  discussed  later  for  realistic 
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geometries)  it  is  convenient  to  model  the  air- 
craft by  a single  geometric  form.  Since  the 
Laplace  equation  is  separable  in  ellipsoidal 
coordinates,  a versatile  3-D  form  that  can  be 
defined  by  one  coordinate  of  a 3-parameter 
ellipsoidal  coordinate  system  is  the  tri-axial 
ellipsoid.  Here  solutions  are  expressed  in 
elliptic  integrals,  that  need  special  tables  or 
computer  programs  for  numerical  evaluation. 
However,  if  we  specialize  the  ellipsoid  to  a bi- 
axial (rotationally  symmetric)  prolate  spheroid 
or  to  an  elliptical  cylinder,  simpler  solutions 
expressed  in  terms  of  trigonometric  functions 
(arc  tangent)  are  obtained.  In  many  cases  these 
simplified  mathematical  models  are  sufficient  to 
study  conditions  on  the  edges  of  the  wings  and 
on  the  tail  (using  the  elliptic  cylinder)  and 
conditions  on  the  fuselage  (using  the  prolate 
spheroid) . When  we  specialize  to  a circular 
cylinder  the  difference  in  the  induced  surface 
fields  due  to  the  separate  orthogonal  components 
of  an  incident  field  vector  (the  cause  of  the 
site  error)  is  still  significant,  as  will  be 
shown.  However,  if  we  specialize  further  to  a 
sphere  this  difference  vanishes  because  of  the 
extreme  symmetry  of  the  sphere  (see  e.g.  Table 
1).  Therefore,  the  sphere  has  no  site  error  and 
cannot  be  used  to  model  the  misdirection  effect 
of  an  airplane  body. 

Consider  an  ellipsoid  representing  an  air- 
plane or  a part  thereof.  Figure  1 shows  the 
ellipsoid  projected  onto  the  x-y  plane  of  a 
cartesian  coordinate  system.  It  is  aligned  with 
the  3 axes,  and  has  semi-axis  lengths  a,  b,  and 
c along  the  x,  y and  z axes,  respectively.  (The 
z axis  is  not  shown.)  We  assume  that  the  top  of 
the  airplane  body  points  in  the  +z  direction, 
and  that  the  incident  magnetic  field  is  parallel 
to  the  x-y  plane.  The  incident  wave  propagates 
with  angle  of  attack  0O  with  respect  to  the 
y-axis  as  shown,  and  with  the  magnetic  field 
lines  perpendicular  to  this  direction.  One  can 
consider  the  magnetic  field  distortions  at  sev- 
eral points,  such  as  A,  B,  C,  and  D in  the  fig- 
ure. These  can  represent  for  example  positions 
along  the  centerline  on  top  of  a fuselage  or 
vertical  stabilizer  and  rudder  ("tail  fin"),  or 
along  a wing  or  horizontal  stabilizer  and  ele- 
vator. 

Later  we  will  consider  the  axially  symmet- 
ric case  in  which  the  y-axis  is  the  axis  of 
rotational  symmetry.  Thus  we  consider  the 
solution  to  be  a superposition  of  the  two  pri- 
mary solutions:  One  is  the  "transverse"  case 

where  0o»O°,  the  Incident  field  being  in  the  x- 
direction,  perpendicular  to  the  axis  of  sym- 
metry. The  other  is  the  "axial"  case  where 
6o»90°,  in  which  the  incident  field  is  in  the 
y-direction,  parallel  to  the  axis  of  symmetry. 

Before  considering  the  spheroid,  let  us 
first  consider  the  limit  in  which  b becomes 
infinite.  The  ellipsoid  then  elongates  into  a 
long  cylinder  or  wire  parallel  to  the  y-axis 
and  of  constant  elliptic  cross-section,  as 
shown  in  Fig.  2. 


CIRCULAR  CYLINDER  (c=a) 

Here  we  treat  the  circular  cylinder,  where 
C“a  and  b=°°,  first  in  the  axial  case  where 
9o“90°.  In  this  case  the  field  is  parallel  to 
the  cylinder  axis,  and  the  solution  becomes 
trivial.  The  field  is  excluded  from  the  interior 
by  a solenoidal  sheet  current  in  the  surface,  in 
the  azimuthal  direction  about  the  axis.  This 
current  produces  an  internal  field  that  cancels 
out  the  incident  field,  but  the  current  produces 
no  external  field.  Hence  there  is  no  distortion 
in  the  axial  case,  9o=90°. 

Next  we  treat  the  circular  cylinder  in  the 


transverse  case,  where  9 


Figure  3a  shows 


the  field  ("flow")  lines  in  this  case.  The 
solution  of  Laplace's  equation  that  satisfies 
the  boundary  conditions  at  infinity  and  on  the 
surface  is  the  potential  function 


B * (r  + a /r)  *cosi () 
o 


where  r denotes  the  cylindrical  radial  coordi- 
nate, a denotes  the  cylinder  radius,  BQ  denotes 
the  magnetic  field  magnitude  at  infinity  (large 
r)  , and  t|>  is  the  azimuthal  angle  about  the  cyl- 
inder axis  (<M)  along  the  direction  of  ffQ) . It 
may  be  readily  verified  that  the  potential 
defined  by  Eq.  (1)  satisfies  Laplace's  equation 
in  cylindrical  coordinates: 


+ Ii£  + J 0 
, 2 r 8r  2 ^,2 
3r  r 3<i 


One  boundary  condition  is  that  as  r becomes 
large , $ approaches 


~ -Bq  r cosy 


fepresenting  the  potential  of  the  constant  field 
0 along  the  x-direction.  On  the  surface  the 
other  boundary  condition  is  that  the  normal  com- 
ponent of  the  gradient  vanishes,  that  is, 


■ +B  (1  - %-)  cosil* 

dt  O 


vanishes  for  all  ♦ when  r>a.  Thus  Eq.  (1) 
satisfies  both  boundary  conditions. 

The  radial  and  azimuthal  field  components, 
Br  and  B^,  are  given  by  the  components  of  the 
gradient , 


" - - B0’(*  ” %)-co.* 

3 x 
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1 _3$ 
r 3 ip 


-B  (1  + -r)  -siniji 
O L 
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whose  planes  are  parallel  to  the  z-axis  and  the 
direction  of  attack  (i.e.,  skewed  with  respect 
to  the  cylinder  axis,  as  shown  in  Fig.  2). 
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The  distorted  field  of  interest  is  obtained  by 
considering  the  x-component  Bx  at  r-a,  namely, 

(B  ) - B cosiji  -B  , sin£ 

x r 
a y 

- Bq  + Bo(2sin2*  - l)(a2/r2) 

- B + (2sin2*  - 1)B  (7) 

o o 

where  the  first  term  on  the  right-hand  side  of 
Eq.  (7)  represents  the  undisturbed  field,  and 
the  second  term  represents  the  distortion 
caused  by  the  cylinder. 

Thus,  the  perturbation  field  (addition  to 
the  incident  field)  goes  from  -B0  on  the  "side" 
looking  into  the  magnetic  field,  to  +BQ  on  the 
"top"  of  the  cylinder  (where  the  dot  in  Fig.  3 
represents  the  probable  position  of  the  crossed 
loop),  as  i (/  increases  from  0 to  ir/2.  Assuming 
the  crossed  loop  to  be  normally  positioned  (at 
the  dot  position)  on  "top",  it  senses  a total 
field  value  2B0,  that  is,  the  incident  value 
enhanced  by  a factor  2. 

Now  we  are  in  a position  to  calculate  the 
misdirection  at  angle  of  attack  0 . On  top  of 
the  cylinder  B^j  ia  enhanced  by  a factor  2, 
while  By  remains  unaffected.  Hence  we  have 


B 

- tan0 


I t4n9o 


(8) 


as  the  definition  of  the  apparent  angle  of  inci- 
dence, 0,  and  therefore 


ELLIPTIC  CYLINDER 

In  the  more  general  elliptic  cylinder  case 
(cy*a) , the  two  principal  fields  are  again  in  the 
x-direction  (0o«O°,  transverse),  and  the  y- 
direction  (8o»90°,  axial).  Again  as  in  the  cir- 
cular cylinder,  the  axial  field  gives  rise  to  no 
perturbation.  In  the  transverse  esse,  however, 
the  field  lines  are  perturbed  as  shown  in  Fig.  3, 
where  in  Figs.  3b  and  3c  the  lines  are  topologi- 
cally similar  to  the  circular  case.  Fig.  3a. 
Figure  3b  (c<&)  can  represent  a flat  wing,  with 
the  crossed  loop  (large  dot  in  the  middle)  sens- 
ing relatively  little  field  distortion.  Figure 
3c  (c>a)  can  represent  a tail  fin,  wnere  the 
crossed  loop  (large  dot  on  the  edge)  censes  a 
large  distortion. 

It  can  be  shown  chat  for  all  cases  in  Fig. 

3 the  misdirection  is  given  by 


MD  = 8 - 0 ” arc  tan  (R  tana  ) - 3 (10) 

o o o 

where  the  "ratio-factor"  R is  defined  by 


R 


a 

a+c 


(ID 


Thus,  R reduces  to  1/2  when  c-a  as  in  Fig.  3a, 
and  Eq.  (10)  reduces  to  Eq.  (9)  for  the  circular 
cylinder.  At  the  middle  of  a flat  wing  (Fig.  3b) 
where  c«a,  R reduces  to  approximately  unity  so 
that  the  distortion  is  minimal.  On  the  edge  of 
a tail  fin  (Fig.  3c)  where  c»a,  R becomes  small 
and  0 becomes  zero  (due  to  the  dominance  of  Bx 
which  becomes  large),  independent  of  0Q.  Thus 
the  distortion  is  equal  to  -0Q  and  becomes  large 
at  angles  near  90°  (except  at  90°) . 
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HD  5 0 - 0 - arc  tan (4  tan 9 ) - 0 (9) 

o 2 o o 

is  the  amount  of  misdirection  or  bearing  error 
at  angle  of  attack  90. 

From  Eq.  (9)  we  deduce  that  the  bearing 
error  is  zero  (i.e.,  8*0O)  when  9Q«0  and  when 
9o«90°.  The  maximum  error  is  -19.5°  (or  +19.5°) 
occurring  when  9Q  is  54.7°  (or  its  supplement 
123.3°).  This  follows  from  the  vanishing  of  the 
derivative  of  MD  with  respect  to  80.  The  varia- 
tion of  bearing  error  with  angle  0Q  given  by 
Eq.  (9)  is  shown  in  the  "needle  limit"  column  of 
Table  1.  These  results  are  consistent  with  ADF 
loop  antenna  calibrations  shown  in  Ref.  (3). 

At  arbitrary  angle  of  attack  9Q,  the  axial 
field  component  gives  rise  to  skin  currents 
flowing  arlmuthally  around  the  cylinder,  while 
the  transverse  field  component  gives  rise  to 
longitudinal  currents,  in  the  + y directions  on 
the  two  x-sldes  of  the  cylinder.  Thus  the 
resultant  skin  current  flow-lines  are  ellipses 


PROLATE  SPHEROID 

A prolate  spheroid  is  a reasonable  model 
for  a fuselage  of  finite  length.  In  this  case 
(referring  to  Fig.  1,  where  y is  now  the  axis  of 
rotation),  the  radii  are  related  by  a-c,  as  in 
the  circular  cylinder,  but  b is  now  finite. 

Hence  the  results  depend  on  the  aspect  ratio 
t-a/b,  which  is  less  than  unity  for  a prolate 
spheroid  (and  greater  than  unity  for  an  oblate 
spheroid) . The  incident  direction  angle  9o”0° 
if  the  direction  of  incidence  is  parallel  to  the 
long  axis,  and  8 "90°  if  the  direction  of  inci- 
dence is  perpendicular  to  the  long  axis. 

It  can  be  shown  that  the  ratio-factor  R in 
Eqs.  (10)  and  (11)  is  now  replaced  by  GR  (see 
Eq.  (IS)  below),  where  G is  defined  by 

22**  2?  22** 

G(y,t)  - (b  -y  ) /(b  -y  + ty)  (12) 


where  t”a/b,  and  y is  the  y-dlstance  of  the 
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surface  point  from  the  mid-point;  and  R may  be 
defined  by 


R(t)  - 0.5  (1  - 2t2  + N)/(l  - N)  (13) 


with  N defined  by 


(14) 


2 

and  s - (1-t  ) . 
given  by 


Thus  the  misdirection  angle  is 


MD  = 0 - 0 » arc  tan  (GR  tan©  ) - 0 (15) 

o o o 


where  G and  R are  given  by  Eqs.  (12)  and  (13). 
Note  that  G denotes  the  cosine  of  the  tilt  angle 
3 of  the  tangent-plane  with  respect  to  the  hori- 
zontal plane.  Note  also  that  G**l  for  y>0  (the 
crossed  loop  is  centered  on  the  fuselage) . 

As  the  aspect  ratio  varies  from  the  "needle" 
limit  t-*0  (an  infinitely-long  circular  cylinder) 
to  the  opposite  limit  t-*-l  (sphere) , the  factor  R 
varies  from  0.5  to  1.0.  Assuming  the  crossed 
loop  to  be  centered  on  the  fuselage  (G-l) , the 
misdirection  angle  MD-A0  is  tabulated  as  a func- 
tion of  0O  in  Table  1,  for  values  of  t*0,  0.1, 
0.2,  0.5  and  1.0.  (The  corresponding  values  of 
R are  0.5,  0.5207,  0.5591,  0.7100  and  1.0.) 

Note  that  MD-0  if  6o”0°  or  90°  in  all  cases,  and 
has  a maximum  magnitude  for  0Q  somewhere  between 
0°  and  90°. 

As  the  aspect  ratio  t increases  the  largest 
MD  magnitude  decreases  from  its  value  of  19.5° 
for  an  infinitely-long  circular  cylinder.  Thus, 
if  the  crossed  loop  is  centered  on  an  F-106 
Delta  Dart,  which  may  be  approximately  modeled 
by  a prolate  spheroid  with  aspect  ratio  t*0.1, 
the  largest  MD  magnitude  is  predicted  to  be 
18.4°  by  Table  1,  and  occurs  at  0O"55°.  If  the 
aircraft  is  a C-130,  characterizable  by  t 
between  0.1  and  0.2  (but  assuming  t-0.2)  the 
largest  MD  is  16.4°,  occurring  at  0-55°.  As  t 
increases  further  (as  the  fuselage  becomes 
thicker  compared  with  its  length) , the  maximum 
MD  decreases  toward  zero,  while  the  correspond- 
ing 0O  moves  slowly  toward  45°.  In  all  cases 
the  errors  Increase  if  the  loop  position  is 
moved  off  center  toward  the  nose  or  tail. 
Reference  (2)  presents  further  tabulations  of 
Eq.  (15)  for  off-center  positions,  such  as  A, 

B,  C,  and  D in  Fig.  1. 

REALISTIC  NUMERICAL  MODELS 

For  detailed  realistic  geometries,  computer 
methods  must  be  employed.  Our  3-D  computer 
model  numerically  solves  the  Laplace  equation, 
in  Integral  form,  by  a method  of  moments,  sub- 
ject to  the  boundary  conditions  of  uniform  field 


at  Infinity,  and  zero  normal  gradient  at  the 
aircraft  surface.  The  aircraft  surface  is 
approximated  by  a large  number  of  small  quadri- 
lateral "panels"  or  "patches",  as  illustrated  in 
Figs.  4 and  5.  (Our  code  was  adapted  from  a 
fluid-flow  code  due  to  J.  L.  Hess.)  It  should 
be  noted  that  the  problems  treated  here  are 
equivalent  to  the  problems  of  er^or  compensation 
encountered  in  the  calibration  of  low-frequency 
ADF  antennas  for  aircraft  (3). 

In  the  light  of  data  obtained  in  1977  and 
1981  by  the  Air  Force  in  flight  tests  of  a 
Stormscope  crossed  loop  (4,  5,  6),  we  applied 
our  computer  code  to  obtain  a preliminary 
assessment  of  the  possible  Influence  of  site 
errors.  The  1977  tests  (4,  5)  involved  a T-39 
Sabreliner,  with  the  instrument  Installed  near 
the  leading  edge  of  (and  on  the  underside  of) 
the  right  wlngtlp.  The  model  portrayed  in  Fig.  4 
shows  the  wing  modeled  reasonably  realistically, 
while  the  fuselage,  whose  detailed  structure 
should  be  unimportant  in  this  case,  is  modeled 
crudely.  The  1981  tests  (6)  involved  a C-130 
Hercules,  with  the  Instrument  located  under  the 
fuselage  near  the  tail  and  close  to  the  cargo 
door.  The  detailed  structure  was  modeled  numer- 
ically as  illustrated  in  Fig.  5.  In  this  case 
we  require  realistic  modeling  of  the  fuselage 
and  tall  structure. 

We  consider  first  the  T-39  results,  and  then 
those  for  the  C-130. 

T-39  MODEL  - The  panels  In  Figs.  4a  and  4b 
are  labelled  by  letters  A-H,  denoting  various 
sections,  with  A and  B on  the  fuselage,  and  C-G 
on  the  wing  and  H on  the  wlngtlp.  Each  section 
has  12  panels,  with  Nos.  1-6  on  the  upper  sur- 
face and  Nos.  7-12  denoting  image  positions  on 
the  under  surface  (with  7 under  6,  8 under  5, 

...,  and  12  under  1).  The  Stormscope  Instrument 
position  is  on  Panel  G-ll,  as  indicated. 

Some  selected  preliminary  results  are  as 
follows,  indicating  panel  location,  maximum  mla- 


direction. 

and  angle  of  incidence  at  which  this 

occurs.  For  each  section  we  give  the  optimum 
location. 

A-9: 

+11°  at  135°  (bottom  of  fuselage 
at  midwing) 

B: 

(no  good  location,  vertical  plane) 

C-3: 

+16°  at  30°  (top  of  wing,  behind 
leading  edge) 

C-10: 

+18°  at  5°  (bottom  of  wing,  behind 
leading  edge) 

D-2: 

+8°  at  15°  (top  of  wing,  adjacent 
to  leading  edge) 

D-ll: 

+12°  at  5°  (bottom  of  wing, 

adjacent  to  leading  edge) 

E-lls 

+5°  at  15°  (bottom  of  wing, 

adjacent  tc  leading  edge) 

E-2- 

+10°  at  160°  (top  of  wing, 

adjacent  to  leading  edge) 

F-5: 

+3°  at  110°  (top  of  wing,  ahead 
of  trailing  edge) 

F-10: 

+3°  at  25°  (bottom  of  wing,  behind 
leading  edge) 
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G-10:  +4°  at  80°  (adjacent  to  wingtip, 

bottom,  behind  leading  edge) 

G-3:  +6°  at  115°  (adjacent  to  wingtip, 

top,  behind  leading  edge) 

G-ll:  -7°  at  170°  (Stormscope  location, 

bottom,  leading  corner  of 
wingtip) 

G-2:  +10°  at  95°  (mirror  of  Stormscope 

location,  top  surface) 

H:  (no  good  location,  vertical  plane) 

The  foregoing  represent  optimum  locations 
(where  the  misdirections  are  minimal) . The  mis- 
directions are  larger  at  other  locations. 

The  following  conclusions  may  be  drawn. 

The  optimum  fuselage  location  is  underneath,  at 
midwing  position.  The  optimum  locations  on  the 
wing  (and  in  fact  on  the  whole  airplane)  are 
near  the  wingtip  and  away  from  the  fuselage, 
either  on  top  and  ahead  of  the  trailing  edge,  or 
underneath  and  behind  the  leading  edge.  The 
actual  Stormscope  location  used  was  a reasonable 
choice  (in  the  absence  of  data  on  site  errors) 
but  could  have  been  improved.  The  reported 
bearing  discrepancies  (5)  are  consistent  with 
the  computed  maximum  misdirection  near  the  wing- 
tip, of  the  order  of  10°. 

It  should  be  mentioned  that  these  figures 
apply  to  an  incident  field  lying  entirely  in  the 
horizontal  plane.  The  presence  of  a vertical 
field  component  would  be  associated  with  a 
slanted  lightning  channel.  In  this  way  the 
error  due  to  slant  would  be  coupled  with  site 
error.  It  if  straightforward  to  include  the 
vertical  component  in  the  calculations,  i.e.  to 
study  the  site  errors  associated  with  slanted 
channels . 

C-130  MODEL  - The  geometry  of  the  C-130  was 
simulated  as  illustrated  in  Fig.  5.  The  loca- 
tion of  the  Stormscope,  just  aft  of  the  cargo 
door,  is  indicated  in  Fig.  5 by  a darkened 
panel.  Some  detailed  small  structures  such  as 
engine  pods  have  been  omitted  in  this  prelimin- 
ary work  on  the  assumption  that  such  small 
structures  are  sufficiently  remote  from  the  par- 
ticular location  of  interest  that  their  contri- 
butions should  be  negligible. 

At  the  Stormscope  location  the  maximum 
bearing  error  or  misdirection  (MD)  is  predicted 
to  be  approximately  20°.  This  error  is  experi- 
enced for  incident  lightning  wave  direction  55° 
(left  or  right)  from  the  forward  direction. 
Although  the  instrument  was  located  on  a flat- 
tened area,  which  would  be  expected  analytically 
(ellipsoid  model)  to  result  in  a smaller  MD,  the 
nearby  tall  structure  has  a strong  influence  on 
the  MD.  With  the  tail  structure  omitted,  the 
calculation  predicts  an  HD  value  of  about  12°. 
Restoring  the  tail  structure  raises  the  MD  to 
20°.  Lower  values  of  the  MD  are  predicted  at 
other  locations  on  the  aircraft  (e.g.  of  order 
4°  Just  ahead  of  the  cargo  door) . The  calcula- 
tions are  still  in  progress  and  further  results 
will  be  reported. 

The  overall  agreement  between  the  analyti- 
cal calculations  (e.g.  using  cylinder  and 
ellipsoid  models)  and  the  numbers  obtained  from 
the  computer  model  is  satisfying.  However,  it 


is  evident  that  more  detailed  information  can  be 
obtained  with  the  computer  model  calculations. 

The  analytical  calculations  are  limited  to  simple 
analytical  shapes,  whereas  with  the  computer 
model  one  can  treat  more  realistic  airplane 
shapes. 

NONVERTICA!  LIGHTNING  CHANNELS 

In  this  section  we  consider  the  error  in 
magnetic  direction  finding  due  to  the  effect  of 
nonvertical  orientation  of  a radiating  lightning 
channel,  modeled  by  a radiating  dipole  tilted 
away  from  the  vertical.  This  effect  has  been 
analyzed  previously  (9,  10),  but  only  with  the 
crossed  loop  detector  at  ground  level.  This 
situation  results  in  important  simplifications. 
The  more  general  case,  where  the  detector  is  not 
on  the  ground  but  airborne,  does  not  appear  to 
have  been  considered  in  the  literature.  We  deal 
here  with  the  more  general  case,  where  both 
source  and  observer  are  above  a conducting 
ground  plane.  The  magnetic  crossed-loop  detector 
is  considered  to  be  otherwise  perfect  (no  per- 
turbations due  to  the  aircraft,  etc.).  The  "mis- 
direction" due  to  induced  skin  currents  was 
treated  in  the  previous  section. 

Referring  to  Fig.  6,  we  consider  the  detec- 
tor to  be  at  altitude  h above  the  ground  plane 
(mounted  on  an  aircraft).  The  lightning  radia- 
tion source  is  represented  by  a dipole  of  vector 
dipole  moment  ?,  located  at  altitude  H above  the 
ground  plane.  Define  a cartesian  coordinate  sys- 
tem, wherein  the  ground  plane  is  represented  by 
the  x-y  plane,  and  the  dipole  is  centered  on  the 
z-axis  at  a vertical  height  z-H  above  the  ground 
plane.  The  detector  is  located  in  the  x-z  plants, 
at  a height  Z"h  above  the  ground  plane,  and  at 
horizontal  distance  x-D  from  the  z-axis  (for 
y**0) . Thus,  the  coordinates  of  the  source  and 
detector  in  this  system  are  (0,  0,  H)  and 
(D,  0,  h),  respectively.  Let  R^  denote  the  dis- 
tance from  the  dipole  to  the  detector,  and  let 
R2  denote  similarly  the  distance  from  the  dipole 
image  (at  depth  H below  the  ground  plane)  to  the 
same  detector.  Then  Rj^  and  are  given  by 
RI-D2+(H-h)2  and  Rj-DZ+(H4h)2.  Let  the  dipole 
be  arbitrarily  oriented,  with  components  P_,  P„, 
and  Pz.  We  can  then  compute  the  magnetic  field 
as  the  superposition  of  contributions  f jom  the 
3 separate  vector  dipoles^of  moments  ?xi,  P J, 
and  P£k,  where  1,  j,  and  k denote  unit  vectors 
in  the  x,  y and  z directions. 

According  to  Hertz'  solution  of  Maxwell's 
equations,  the  electric  and  magnetic  fields  due 
to  a radiating  dipole  may  be  expressed  in  terms 
of  a dipole  moment  vector  function  (11,  12).  If 
we  denote  the  Hertzian  time-dependent  function 
by  ?(t-R/c),  the  vector  potential  I may  be  writ- 
ten as 

t - i t (t-R/c)  (16) 

where  the  argument  denotes  the  retarded  time  at 
distance  R from  the  dipole,  and  the  dot  denotes 
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differentiation  with  respect  to  the  time.  Here, 
the  units  adopted  are  Gaussian  and  ? is  expressed 
in  emu,  i.e.,  abamp-cm,  with  R in  cm.  (In  MKS 
units,  an  additional  factor  pQ/4 ir  appears  on  the 
right  side,  where  is  the  magnetic  permeability 
of  free  space,  AmclO  ^ henry/meter;  then  ? is 
expressed  in  ampere-meters,  with  R in  meters.) 

Now  the  magnetic  field  intensity  B is  given 
by 

if  - curl  £(R)  - (R/R)  x (dX/dR)  (17) 

since  A is  a function  of  R only.  Differentiation 
yields 
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(22) 


where  we  identify  the  unsubscripted  variables 
with  the  original  dipole.  It  follows  from  Eqs. 
(19)  and  (21)  that,  at  the  observer's  position, 
the  magnetic  field  of  the  original  dipole  is 
given  by 


“ f-<H-h)Mly]  i 

+ [(H-hJM^  + D Mlz]  j (23) 


dX  t L 

dR  “ r2  Rc 


(18) 


where  the  second  term  results  from  the  dependence 
of  the  time-like  argument  on  R,  and  we  obtain 


+ [-D  Kly]  * 

while  the  magnetic  field  due  to  the  image  is 
given  by 

t2  - t(H+h)M2y]  i 


B - M(R)  x R 


(19) 


+ [-(H+h)^  + Dy  j (24) 


with 


+ [-D  M^]  k 


(20) 


so  that  a = 1/R^  and  B = 1/cR^. 

In  terms  of  our  cartesian  coordinates,  let 


Hence,  the  sum  of  the  two  fields  is: 


if  - if. 


“ t(M2y  - Mly)H  + (M2y  + Mly)h]  1 


“l  ■ Mlxi  + V + ”1*" 

^2  ■ M2xi  + V + *2*' 

R - Di  - (H-h)k 


(21) 


+ [-(M2x  - Mlx)H  - (M2x  + M^h  (25) 
+ (M2x  + Mlz)D]  j 
+ [-(M2y  + Mly)D]  k 


Rj  - Di  + (H+h)k 

where  1,  j and  k denote  unit  vectors  in  the  x,  y 
and  e directions;  and  and  i?2  denote  M^(R^)  of 
the  original  dipole,  and  M2(R2)  refers  to  its 
image,  respectively.  In  this  case,  the  compo- 
nents of  t2>  ^2  are  related  t0  those  ^ 

by: 


For  a vertical  dipole,  all  components  of 
and  i?2  vanish  except  for  M^z  and  M2z.  Then  J 
has  only  a y component,  and  there  is  no  misdirec- 
tion. When  f has  an  x-component  also,  the  mis- 
direction angle  is  given  by  the  src  tangent  of 

-B  /B  (Fig.  6b).  Consider  the  ratio  -B  /B  , 
x y ° x y 

obtained  from  Eq.  (25): 


K -<h2.-mi.)h  - 

*, ' W - «,h  + 
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Using  Eqs.  (20)  and  (22),  we  may  rewrite  Eq.  (26) 


B p H + q h 

_ _*  - 1 L JL 

B pH+qh+pD 
y rx  ^x 


where 


(27) 


"S 


Q 5 


■r:*n  i w.wnf *■!».«!  « 


(<»2  “ a^P  + (B2  ~ Bj^P 
• * • 
(aj  + aj)P  +(62+  BjJP 


(31) 


(1/R23  - i/r13)p  + (i/r22  - i/r12)p/c 
(1/R23  -*  /R^P  + (1/R22  + l/R^P/c 


Px  - (o2  + o1)Px  + (B2  + 61)PX 


py  * <«2  + “l^y  + (62  + 61)PV 


px  - (a2  + 0l)Px  + (e2  + B1)PZ  (28) 


qx  - (a2  - 0l)Px  + (B2  - 81)PX 


qy  " (a2  " °l)Py  + (B2  “ Bl)Py 

• ••  • ••  » •• 

and  where  P , P , P , P , P , P refer  to  the 

x x y’  y 1 t 

original  dipole  components.  Divide  through 
numerator  and  denominator  by  p^D,  and  note  that 
we  may  write: 


In  the  case  of  harmonic  radiation,  of  angu- 
lar frequency  w,  we  may  replace  P by  juP  in  Eq. 
(31),  so  that  Q is  given  by  the  conr,  .ex  quantity, 
with  k=u/c: 


, 3 3.  ....  2 2, 

(*1  - x2)  + jk(x1  - x2) 


(x3  + x3)  + jk(x2  + x2) 


(32) 


«...  ± JL. 

a'  + jt/ 


where 
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are  the  original  dipole  moment 
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vector  components,  and  where  P and  P are  the 
time  derivatives  of  the  vector  magnitude  F. 
then  Px/Px.  Py/pz*  1X/P2»  *nd  1y/P2  become, 
using  Eqs.  (28)  and  (29): 


px  Px 

»/P, 


-S  - Q-* 
P.  P_ 


P P 
JL  m JL 
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0 JL 
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■ Re  Q + j Im  Q 


Finally,  the  ratio  -B  /B  may  be  written: 

* y 


.i 


(33) 


Bx  Ps 


s tan  (MD) 


(34) 


(30) 


y 1 + Gp 


where  HD  is  the  misdirection  angle  and  where 
G is  a geometrical  factor: 


where 


a H + Qh 
G - D 


(35) 
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If  the  detector  is  on  the  ground  plane 
(h«0) , the  problem  simplifies  greatly,  and  we 
have  a^"02  and  3^*$2>  80  that  Q“0.  Then  G 
reduces  simply  to  G”H/D,  and  we  have 


1 + h!x 
DP 

z 


which  is  identical  to  the  result  given  by  Uman 
et  al  (Ref.  10  - see  plots),  re-expressed  in  our 
terms,  for  a detector  on  the  ground  plane.  An 
equivalent  result  was  also  given  earlier  for  this 
case  by  Kalakowsky  and  Lewis  (9) . 

If  the  ground  plane  is  completely  noncon- 
ducting (no  image),  then  the  case  of  the  detector 
at  altitude  h above  ground  also  simplifies 
greatly.  In  this  case  a and  8 vanish,  and 
Q—l.  Then  1 1 


D P 


, . (H-h)  X 
DP 

z 


That  is,  H is  replaced  in  Eq.  (36)  by  the  "net 
altitude,"  H-h. 

A geometric  interpretation  of  Eq.  (36)  (or 
(37))  is  obtained  by  rewriting  Eq.  (34)  in  the 
form,  assuming  Q— 0 : 


P P 

(D  + H ~)  B + H =*  B 
' P x P j 

z z 


This  expression  implies  that  the  "misdirection," 
in  the  plane  of  the  observer  and  perpendicular 
to  S at  the  observer,  passes  through  the  point 
whose  x,  y,  z coordinates  are  -(P„/P,)H, 
-(Py/PS)H,  and  zero.  This  point  (indicated  by 
the7 letter  I in  Fig.  6a)  is  alto  the  intercept 
of  the'  extended  line  of  the  dipole  with  the  z»0 
plane.  This  simple  geometric  interpretation  for 
a detector  on  a conducting  ground  plane  was 
pointed  out  by  Kalakowsky  and  Levis  (9). 

The  correction  factor  for  finite  aircraft 
altitude  is  given  by  the  ten#  Qh  in  Eq.  (35). 
Assuming  a frequency  of  50  kHz  (Stormacope  cen- 
tral frequency),  a lightning  distance  15  km 
(■D),  a lightning  altitude  3 km  (»H) , and  air- 
craft altitude  3 km  (»h) , we  obtain  from  Eq. 

(33)  the  approximate  value  -1/3  for  Re  Q.  Thus, 
for  this  case  the  geometrical  factor  G is 
reduced  to  about  2/15  compared  with  its  zero- 
altitude  (h«0)  value  of  3/15. 


APPENDIX.  USE  OF  ELECTRIC  ANTENNA  TO  RESOLVE 
180-DEGREE  AMBIGUITY 

A simple  way  of  seeing  how  correlation  of 
the  phases  of  the  electric  and  magnetic  vectors 
can  be  used  to  resolve  the  ambiguity  is  the  fol- 
lowing. 

Consider  one  loop  of  a narrowband  crossed- 
loop  system  aligned  with  two  sources,  one  in 
"front"  and  one  in  "back".  In  Fig.  A-l,  the  plane 
of  the  loop  is  perpendicular  to  the  plane  of  the 
paper  and  is  aligned  with  the  sources  (upper  dot 
■ front,  lower  dot  ■ back).  The  sources  are 
assumed  to  be  electric  dipoles  in  the  same  plane 
as  the  detector,  but  oriented  perpendicular  to 
the  plane  of  the  paper.  Four  sets  of  vectors 
are  shown,  at  (a),  (b) , (c)  and  (d) , with  E,  H, 
and  P in  each  set  denoting,  respectively,  the 
electric  vector,  magnetic  vector,  and  direction- 
of-propagation  vector,  respectively.  Sets  (a) 
and  (b)  correspond  to  radiation  from  the  front, 
while  sets  (c)  and  (d)  correspond  to  radiation 
from  the  back.  Let  signs  be  associated  with  E 
and  H,  such  that  E pointing  upward  and  downward 
denotes  a positive  and  negative  amplitude, 
respectively.  While  H pointing  right  and  left 
denotes  also  a positive  and  negative  amplitude, 
respectively  (readily  convertible  into  loop  volt- 
age signs).  The  four  possibilities  may  be  tabu- 
lated as  follows: 

E H Source 


It  is  evident  from  the  table  that  E and  H have 
opposite  signs  if  the  source  is  in  front,  while 
they  have  the  same  signs  if  the  source  is  in 
back.  It  is  not  difficult  to  apply  the  same 
ptweedure  to  the  second  loop. 

Thus  it  is  readily  seen  why  electrical  noise 
can  produce  180°  errors.  If  the  signal- to-noise 
ratio  is  too  low  or  the  electronics  is  ineffec- 
tive, electrical  noise  pulses  of  the  wrong  sign 
can  overwhelm  the  electric  field  amplitude  and 
make  it  appear,  say,  negative  when  it  should  be 
positive,  and  vice  versa.  Similar  errors  may  be 
caused  by  unintentional  phase  shifts  in  the  sig- 
nal amplifiers. 
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Fig.  A-l  - Resolution  of  the  180-degree  ambigu- 
ity 


Fig.  3 - Transverse  field  lines  around  circular 
end  elliptic  cylinders  (large  dot  indicates 
detector  position) 
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